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Abstract: Applying the AdS/CFT correspondence, the expansion of the heavy-quark potential of 
N — A supersymmetric Yang-Mills theory at large Nc is carried out to the sub-leading term in the 
large 't Hooft coupling at zero temperature. The strong coupling corresponds to the semi-classical 
expansion of the string-sigma model, the gravity dual of the Wilson loop operator, with the sub- 
leading term expressed in terms of functional determinants of fluctuations. The singularities of 
these determinants are examined and their contributions are evaluated numerically. 
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1. Introduction 

AdS/CFT duality ||, | remains an active field of research. Motivated by the isomorphism 
between the isometry group of AdSs and the conformal group in four dimensions, it was conjectured 
by Maldacena that a string theory in AdSs x corresponds to a four dimensional conformal field 
theory on the boundary. A prominent implication of the conjecture is the correspondence between 
the type JIB superstring theory formulated on AdSs x and A/" = 4 supersymmetric Yang-Mills 
theory (SYM) with the isometry group 0(6) of dual to the R-symmetriy group SU{A) of SYM. 
In particular, the supergravity limit of the string theory corresponds to the leading behavior of 
SYM at large Nc and large 't Hooft coupling 

A = 5ym^^c = ^. (1.1) 

with L the AdS radius and a' the reciprocal of the string tension. This relation thereby opens 
a new avenue to explore the strong coupling properties of SYM and sheds new lights on strongly 
coupled QGP created in RHIC in spite of the difference between SYM and QCD. Among notable 
successes on the RHIC phenomenology are the equation of state the viscosity ratio |^ and jet 
quenching parameters as well as the energy lossQ. 

The heavy quark potential (the potential energy between a heavy quark and its ant i- particle) 
of QCD is an important quantity that probes the confinement mechanism in the hadronic phase 
and the meson melting in the plasma phase. It is extracted from the expectation of a Wilson loop 
operator, which can be measured on a lattice. In the case of A/" = 4 SYM, the AdS/CFT duality 
relates the Wilson loop expectation value to the path integral of the string-sigma action developed 
in Ref for the worldshect in the AdSs x bulk spanned by the loop on the boundary. To 
the leading order of strong coupling, the path integral is given by its classical limit, which is the 



minimum area of the world sheet. From the Wilson loop of a pair of parallel lines, Maldacena 
extracted the potential fmrction in J\f = i SYM at zero temperature jl^, 

nr)=-^^^ -0.2285^ {1.2) 

with r the distance between the quark and the antiquark. Introducing a black hole in AdS bulk, 
the potential at nonzero temperature as well as that for moving quarks have been obtained by a 
number of authors . The field theoretic aspects of the potential (L2) and its finite temper- 



ature counterpart as well as their implications on RHIC physics were discussed in Ref.|13] |14| |12l 



As was pointed out in Ref.llQ], the "heavy quarks" underlying the Wilson loop (1.2) in J\f = 4 



SYM are actually heavy W bosons resulted in a Higgs mechanism, which implement the funda- 



mental representation of SU{Nc). Since the function (1.2) measures the force between two static 
fundamental color objects, we shall borrow the terminology of QCD by naming it the heavy quark 
potential throughout this paper. 

The strong coupling expansion of the SYM Wilson loop corresponds to the semi-classical ex- 
pansion of the string-sigma action and reads 

n.)^-^^[i.;^.o(i)] 

for the heavy quark potential. Computing the coefficient k is the main subject of the present paper. 
AC comes from the one loop effective action of the world sheet fiuctuations around its minimum area. 
This effective action has been obtained explicitly for some simple Wilson loops including parallel 



lines 15 [16 and is expressed in terms of functional determinants. Evaluating these determinants. 



we end up with the numerical value of k, 

K~ -1.33460. (1.4) 

The classical solution of the string-sigma model and the one loop effective action underlying k 
is briefly reviewed in the next section. There we also outline our strategy of computation, which 
is along the line suggested in [|l6|. We parametrize the string world sheet of the single Wilson 
line or parallel lines by conformal coordinates. Then a scaling transformation is made that leaves 
the measure of the spectral problem of the functional determinants trivial. Instead of solving the 
eigenvalue problem of the operators underlying the determinants, we use the method employed in 
]l7t , which amounts to solve a set of ordinary differential equations. Unlike the straight Wilson line 
and the circular Wilson loop dealt with in |l^ , some of differential equations for the parallel lines 
are not analytically tractable. The presence of various singularities makes numerical works highly 
nontrivial. It is critical to isolate the singularities analytically in order to obtain a robust numerical 
result. So we did and the procedure is described in sections 3 and 4. The finite terms of the scaling 
transformation of the determinants involved are examined in section 5 and we find them adding up 
to zero. In section 6, we discuss our results along with few open questions. Some technical details 
are explained in appendices. Throughout the paper, we shall work with Euclidean signature with 
the AdS radius L set to one. 

2. The one-loop effective action 



Let us begin with a brief review of the classical limit that leads to the leading order potential (|l.2| ) 
The string-sigma action in this limit reduces to the Nambu-Goto action 
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with g the determinant of the induced metric on the string world sheet embedded in the target 
space, i.e. 



9aP 



= G 



(2.2) 



where and G^^ are the target space coordinates and the metric, and cr" with (a — 0, 1) 
parametrize the world sheet. The target space here is AdSs x , whose metric may be written as 



(2.3) 



with dr^s the element of the solid angle of S^. The physical 3-brane resides on the AdS boundary 
z = 0. The string world sheets considered in this paper are all projected onto a point of in the 
classical limit. 

The Wilson loop of a static heavy quark, denoted by Ci, is a straight line winding up the 
Euclidean time periodically at the AdS boundary. The corresponding world sheet in the AdS bulk 
can be parametrized by t and z with x constant and extends all the way to AdS horizon, z — > oo. 
The induced metric is that of AdS2, given by 



ds^\Ci] 



dz^) 



with the scalar curvature 



R 



-2. 



Substituting the metric ( |2.4| ) into (|2.l| ), we find the self-energy of the heavy quark 



dz 

72- 



(2.4) 
(2.5) 

(2.6) 



with T — )• (X) the time period. Notice that we have pulled the physical brane slightly off the 
boundary to the radial coordinate z = (5, as a regularization of the divergence pertaining the lower 
limit of the integral (2.6). 

The total energy of a pair of a heavy quark and a heavy antiquark separated by a distance r, 
can be extracted from the Wilson loop consisting of two parallel lines each winding up the Euclidean 
time at the boundary. This Wilson loop will be denoted by C2 and the world sheet in the bulk can 
be parametrized by t and z with x^ = ^(z) and x'^,x^ — const.. The function ^(z) is determined by 
substituting the induced metric 



ds'[C2] = -^{dt^ 



1 



dz"^}, 



into the action (2.1) and minimizing it. We have 



dz'- 



V4 



(2.7) 



(2.8) 



The maximum bulk extension of the world sheet, zq, is determined by the distance r between the 
two lines at the boundary and we find that 



£!(i), 

(2^)i 



Substituting (2.8) into (2.7), we end up with the induced metric 



ds^ps] = — [dt 



rdz' 



(2.9) 



(2.10) 
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and the scalar curvature 



i? = -2 

The energy of the heavy quark pair is therefore given by, 

1 1 Z"^" 

E[C2] = 7^5ng[C2] = —zl 



dz 



(2.11) 



(2.12) 



where the same regularization is apphed to the lower limit of the integral. 

The heavy quark potential is obtained by subtracting from ( 2.12| ) the self energy of each 
quark(antiquark) , i.e. 



V 



\im{E[C,]-2E[Ci]) = — 
s^o+ Tra' 



1 



(2.13) 



and is divergence free. Carrying out the integral and substituting in the relations (2.9), we derive 

dH). 

The one loop effective action, W is obtained by expanding the string-sigma action of Ref . Q to 
the quadratic order of the fluctuating coordinates around the minimum area and carrying out the 
path integral|15 [16|. We have 

W[Ci] = - In 

for the static quark or antiquark and 

W[C2] = - In 



det^(-i7"Va +T3) 



det^(-V2 + 2)det5(-V2) 



det^(-i7"Vc«+T3) 



(2.14) 



det5 (-V2 + 4 + i?) det(-V2 + 2) det^ (-V^) 



(2.15) 



for the quark pair. The determinants in the denominators of (2.14) and (2.15) come from the 
fluctuations of three transverse coordinates of the AdS sector and five coordinates of with the 



Laplacian given by the metric (2.4) or (2.10). The determinants in the numerators come from the 
fermionic fluctuations, where we have introduced 2d gamma matrices, 7o = 7'^ = ''"2j 7i = 7^ = ti 
and 7o7i = —iT^ with ti, T2 and T3 the three Pauli matrices. In terms of the zweibein of the world 
sheet, e^, we have 7^ = SaJj with j — 0,1 and the covariant derivative 



d 1 



(2.16) 



9cr" ' 8' 

with ut'^J the spin connection corresponding to ( 2^ ) or ( |2.10| ) . The power " 4" comes from eight 2d 
Majorana fermions each of which contributes a power 1/2. The one loop correction to the heavy 
quark potential is then 

AV^ lim ^ lim (W^[C2] - 2W^[Ci]). (2.17) 

The effective action VF[Ci] or W[C2] suffers from the usual logarithmic UV divergence, which is 
proportional to the volume part of the Euler character 



dtdzy/gR 



(2.18) 



of each world sheet with the same coefficient of proportionality ||l q| . It follows from ( [2.4| ), (2.5), 
(|2.10|) and ( 2.11 ) that the integral ( 2.18 ) for the parallel lines is exactly twice of that for the single 
line in the limit 6^0. We have indeed that 



J d^a^R^T 



dz , , 2T 



(2.19) 
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for the single line and 



(fay/gR = 2T J dz 



AT 



(2.20) 



for the parallel lines. Therefore the UV divergence as well as the conformal anomaly cancel in 
the combination of ( 2.17 ) in the limit 5^0. As a contrast, the volume integral J (Pa^/g of the 
parallel lines differs from twice of that of a straight line by a finite quantity in the same limit. The 
UV divergence associated to the volume integral cancels within each effective action of ( ^.14 ) and 
(2.15)[|8j. Furthermore the limit 5 — >■ 0+ of the UV finite term of ( [2.17| ) also exists as we shall see. 
The world sheet of the parallel lines covers the coordinate patch {t, z) twice, which gives rise to 



an artificial singularity of the Laplacian's in (2.15) at z — zq and adds difficulties to the numerical 
works. To avoid the problem, we shall work with a conformal coordinate patch (r, cr) that the world 
sheet (2.10) covers only once. This is also suggested in 16 . The new coordinates involve Jacobi 
elliptic functions Jl9[ pO[|of modulo fc = -1= and are defined by 



z = ZQCncr t ~ —-=t 
v2 



In terms of the new coordinates, the metric ( 2.1C ) takes the form 

1 



2cn2(7 



(dr^ + da-" 



and the scalar curvature ( p.ll ) becomes 



R= -2(1 + en-* cr) 



(2.21) 



(2.22) 



(2.23) 



The nonzero component of the spin connection with cartesian indexes (0,1) referring to the coordi- 
nate differentials dr and da reads 



,01 



sncrdncT 



(2.24) 



We shall use the the same time variable r to describe the world sheet of the straight line and rescale 



the z coordinate by z = -^C: leaving the conformal structure of (2^) intact, i.e. 

ds^[Ci] ^ ^{dr^ + dC^). 



(2.25) 



The spin connection corresponding to ( 2.24 ) is given by uj^ 



,01 



,10 



The range of each 

coordinate variable is — ^ < r < ^, —K < a < K and < C < oo where T = —T and K is the 
complete elliptic integral of the first kind. 



r2 (^) 
K = ^ ~ 1.8541. 



The operators underlying the determinants of (2.14) are given explicitly by 



Ao[Ci] 
Ai[Ci] = - 



C'Ao[Ci 



2 = -C(^ + ^)+2^C^Ai[Ci 



(2.26) 

(2.27) 
(2.28) 
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and 



Df[Ci] = -z7"Vc« + T3 - -«C hj7 - ^ ^1 - <^^2 + T3 = C^f[Ci] 



d 



'dr 



(2.29) 



Similarly, the explicit expressions of the operators underlying the determinants of ( p.ll ) reads 

-V^ = -2cn2CT 



Ai[C2] = - 



-2cn^cr 



& 



A2P2] = -V2 + 4 + i? = -2cn2c7 + + 2(1 - cnV) = 2cn2aA2[C2] 



dr 



^)^2cnVAo[C2], 
— j + 2 = 2cn2CTAi[C2], 



and 

Df[C2] 



iV2c 



d_ 



sncrdncrX 
2cn(T / 



Ti - i\/2cn(T— r2 + ra ee cii(7Df[C2] 

OT 



(2.30) 

(2.31) 
(2.32) 

(2.33) 



The difference between the operators with hats and those without hats is the measure of the spectral 
problem defined by them. While the measure is trivial with respect to the operators with hats, 
changing the measure may introduce additional terms to the logarithm of each determinant and 
their contribution will be examined in section V. For this reason, the effective action is decomposed 
into two pieces, i.e. W[Ci] = W'i[Ci] + 1^2[Ci] for the single Wilson line and 1^[C2] = VFi [C2] + ^^2^2] 
for the parallel lines. We define 



Wi[Ci 
W2[Ci] = -4 In 

Wi[C2] = - 



= —^1^ 5"^^ 3-^^ 

det^ Ao[Ci] det^ Ai[Ci] 

detDF[Ci]\ 2 Ao[Ci] 
det'^ Df[C2] 



^ In Mil 
2 ""AiPi] 



In 



and 



det5 A0P2] dot A1P2] dct2 A2[C2 
Ai[C2 



M.2[C2]^-41n^i^ia« + ^ln^°[^^ 
\detDF[C2]\ 2 



In- 



A0P2] Ai[C2] 
Correspondingly, the coefficient k defined in (^^) is given by k = ki 



1, A2P2 
- In — 

2 A2[C2 



K2 with 



and 



Ki 



K2 



lim {Wi[C2]-2Wi[Ci]), 



lim (W^2[C2]-2VK2[Ci]), 



(2.34) 

(2.35) 
(2.36) 

(2.37) 

(2.38) 
(2.39) 



where we have used the relation between T and T and converted zq to r via (2.9). 

Making a Fourier transformation of the time variable r, each functional determinant of (2.14) 
and ( 2.15| ) is factorized as an infinite product of its Fourier components with each Fourier component 
obtained by replacing the time derivative ^ in A's of (2.27)-(2.33) by — with uj a frequency 
variable. Substituting the Fourier product of ( 2.14 ) and that of ( 2.15 ) into ( 2.38 ), we find that 



Kl 



°° duj_ 7^2(cJ) 

_^ 2^ nliuj) 



(2.40) 
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The functions Tli{uj) and Tl2{io) are the Fourier components of the determinant ratios of ( 2.14 ) 

det2?^(tj)det2?2 (w) 



and ( 2.15 ) , given by 



and 



7^2 (w) 



detV§ (a;)detX'f (w) 
deti:>2 (a;)dctD2 (-w) 



detDg" (a;)detDi(w)deti:'| (w) 
where the Fourier transformation of the operators A's and Dp's are given by 



(2.41) 



(2.42) 



dC 



da^ 



d^ 

d<T^ 



cn^cr, 



cn-'CT 



and 



1 =p v^sncrdncr \ 



(2.43) 
(2.44) 
(2.45) 
(2.46) 
(2.47) 

(2.48) 



2cn2(T j 

Let us explain the transformation we made on the fermionic determinants det [Ci] and det Dp [C2] 
Replacing the time derivatives in (2.29) and (2.33) by — ioj, we find that 



Df[Ci] = -i\-rT--^]Ti- UJT2 + -T3 



1 



and 



Df[C2] = -i\ — 
It is straightforward to verify that 



dC 2C 



d sncrdncr \ 
da 2cn(T J 



c 



Tl - LJT2 



V2( 
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cncr 



and 



DUCi] = yCC/diag.(I?+(c.),2?_(c.))[/ti= 



Dl[C2] = V^Udiag.{D+{iu),D^{uj))U^^=, 



(2.49) 



(2.50) 



(2.51) 



(2.52) 



where ?7 is a 2x 2 matrix that diagonalizes T2, 'D±{uj) and D±{uj) are given above by (2.43) (2.44) and 
( |2.48D . Therefore dei'^ DplCi] = det^ r'+(w) det^ 2?_(w) and Aei'^ bF[C2] = det^ i:)+(w) det^ i:)_(a;). 
The evaluation of the integral ( 2.4Cl|) will be discussed in the next two sections. 
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3. The evaluation of ki 



— - analytical part 



Evaluating a functional determinant stemming from a semi-classical approximation to a quantum 
mechanical system is an old subject of many research works. In one dimension a short cut was 
discovered by a number of authors p2[ that does not require a solution to the spectrum problem 
involved. Consider two functional operators 



Ha 



+ Va{x) 



(3.1) 



with a = 1,2 defined in the domain a < x < b where Va{x) > under the Dirichlet boundary 
condition, it was shown that the determinant ratio 



dcti/2 f2{b\a) 



(3.2) 



detHi hib\a) 
where faix\a) is the solution of the homogeneous equation 

Ha(f> = 0, (3.3) 
and /4(0|a) — 1. In terms of a pair of linearly independent 

Ua{a)Va{b) - Ua{b)Va{a) 



subject to the conditions fa 

solutions of (|3.3D, {Ua,Va), 



Ub\a) 



(3.4) 



W[Ua,Va] 

where the Wronskian is x-independent. With appropriate modification of the conditions 

imposed on fa{x\a), the formula (3.2) can be tailored to cover other boundary conditions. This 
method has been employed recently in |17 to calculate the one loop effective action of the single 
line Ci or that of a circular Wilson loop. See |2^ for a review on other applications. 

Coming back to the semi-classical correction of the heavy quark potential, the operator Ha 
corresponds to one of the operators ( 2.43 ) -( 2.48 ). We shall retain {u,v) for a pair of linearly 
independent solutions of the homogeneous equation ( [3.3| ) with Ha given by an operator pertaining 
to the single line and denote that of the corresponding equation of the parallel lines by (j],^). Eq. 
(3.3) with Ha given by an operator of ( 2.43 )-( 2!45| ) can be solved analytically and we may choose 
the following pairs of independent solutions 



ui — coshw^ — 



uq = sinhwC = uo{ujQ 
sinhcjC 



= ui(wC) 
u+ = ui{ujC) 



Wo = 
vi=\l 



1 



wo(wC), 



v+ = vi{u}C) 



and 



r]o = sinh[j(T = coshcjCT. 

with their Wronskian's all given by 

W[uo,vo] = VK[r/o,Co] = W[ui,vi] = W[u±,v±] 

The equations 1^ with Ha given by ( pl6| ), ( pl7| ) and (|1|), 



= 0, 



(3.5) 

(3.6) 

(3.7) 
(3.8) 

(3.9) 
(3.10) 

(3.11) 
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and 



cn^CT 



0, 



1 =p v^sncrdncr \ 
2cn2cr ) 



(3.12) 



(3.13) 



do not admit analytical solutions for a; 7^ 0. Eqs.(3.11) and ( ^.12 ) have a = ±K as regular points 
with the same pair of indexes (2,-1) there. The equation D±{uj)(j) = has a regular point a = 
with the indexes (2,-1) and a = ±K is an ordinary point of it. We associate 77's to the vanishing 
solution at (T = —K and ^'s to the vanishing solution at a — K with the normalization conditions 



lim 



and 



lim 



lim 



lim 



lim 



6(.) _ 



(3.14) 



1. 



a^K uj^{K -a)'^ a^Kuj^{K-aY a^Koj'^{K-aY 
Furthermore, we require 

i^-{-K)=^+{K)=Q 

and 

Tl'_{-K)=.~i'^{K)=uj 

with the prime the derivative with respect a. On account of the eveness of -Di(a;) and D2{uj) with 
respect to cr, we have 

6,2(^7) =?7i,2(-a). (3.18) 
It follows from the relation between D+{uj) and D^{uj) that 



(3.15) 

(3.16) 
(3.17) 



r,_(r7)=e+(-a) C^{a) = r,+ i-a) 



(3.19) 



Each differential equation of (3.11), (3.12) and (3.13) is of the form of an one dimensional Schroedinger 
equation in a non negative potential at zero energy and does not admit a bound state subject to 
the Dirichlet boundary condition. Therefore we expect that 



as a ^ K and 



^71,2(0-) 
ei,2(a) 



_ Cl,2(^) 

' uj{K-a) 

ui{K - a) 

_ Claim) 
' uj{K + (j) 



(3.20) 
(3.21) 

(3.22) 
(3.23) 



as (T — !• —K. The coefficients of divergence, Ci{uj), C2(w) and C±{oj) are related to the Wronskian's 



via 



with j — 1, 2, ±. We have, in addition, 

V+{K) = 



3w 



U-K) 



(3.24) 



(3.25) 
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It follows from the symmetry property ( ^.19 ) that 

C+ (cj) = C_ (w) = C3 (w) . (3.26) 

For uj >> I, the solutions rj's and ^'s can be approximated by WKB method and we find the 
asymptotic forms 

Ci (a.) = ^e2^"(l--^ + ...), (3.27) 

and 



C2(c.) = |e^«-(1- - + ...), (3.28) 



where the constant 



C3(u;) = V^-(l-^ + ...), (3.29) 
2ttI 



0.84721 (3.30) 



The details of the derivation are shown in the appendix A. The small uj behavior can be obtained 
by introducing an alternative set of solutions, normalized differently, 

m,2,+ = , ^i2-cr) = 2 ^-3^ 

and 

77_(cr) = C+(o-) = • (3-32) 

Defining the coefficients C"s by the diverging behavior 

m.2,-{<y) ^ —f7 + -, (3.33) 

as a ^ K and 

c i \ ^l,2,+ (^) , /Q Q/|^ 

fl,2,+ cr) — — + ..., (3.34) 

A + (J 

ascr^ -J-r,wefindthatCi,2(w) = ^3(71,2(0;) and C3(w) = C±(a;) = a;2C'±(w). Since (71,2(0), Ci (0) 7^ 
and are well defined (determined by eqs.( |3.11 )-(3.13) at w = 0, see Appendix B for details) we 
have the small w behavior, 

247r^ , , 

Ci(w) ~ — rryw^ ~ 10.166557w^ (3.35) 

C2(w) ~ ^TTT'^ - 5.0832785w^ (3.36) 



and 



C3(u;)~4c^2^ (3.37) 



In the regularized version, the physical brane, located at z = S cut the world sheet of the 
parallel lines at —K + e and K — e with 

5 = ZQcn{K - e) ~ ^e. (3.38) 
v2 

In another word, the domain of a coordinate is [—K + e, K — e] under the regularization, and we 
shall impose the Dirichlet boundary condition there. The corresponding domain of the single line 
becomes e < ( < Z with Z a large ( cutoff which will be set to infinity at the end. Designate 
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Uaioj) to the quantity (3.4) of the single Wilson line and Uaii^) to that of the parallel lines with 
a = 0, 1, 2, ± corresponding to the indexes of the operators ( 2.43| )-( 2.48 ), we have ||2^ 



7^l(w) 



Ul{u)Ui{uj) 
uI{lo)VJ{^) 



and 



7^2(cJ) = 



Ul{u)Ul{uj) 



The last step of ( |3.39 ) follows from the solutions ( ^.5[ )-(p^), which imply that 



2uj 



and 



2a; 



as Z — cxD. It follows from ([3.g|) that 



Uoioj) = 



sinh2(i^- e)o 



(3.39) 
(3.40) 

(3.41) 
(3.42) 

(3.43) 



The symmetry ( 3.19 ) implies that U+{uj) = U-{uj). 

To proceed, let us introduce wq that satisfies the inequality 1 << << j and divide the 
integral ( [2.4^ ) into two terms, with 



and 



For the integrand of (3.44), we may approximate 



duj In 



1 



and 



U±{u:) 



lUJt 

3a;3e2 



(3.44) 
(3.45) 

(3.46) 
(3.47) 
(3.48) 



Only one term of the numerator of (3.4) contributes to each case of (3.47) and ( ^.48 ) and the other 
term is suppressed by a power of e. Together with ( 3.43| ), we obtain that 



with 



duj In p{uj) 



(3.49) 



(3.50) 



Ci(cj)C|(w) sinh^ 2Ku; 

and the approximation becomes exact in the limit e 0. It follows from the asymptotic behaviors 
(|I27|) , (|^ and (|^ that 

\np{uj) = o(^j (3.51) 
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for cj >> 1 and 



. , 64\/2 

~ —fTTTu: ~ 2.19171a. 



(3.52) 



as w — >■ 0. The very fact that the integration of k< converges in the limit — > oo indicate that k> 
vanishes under the same hmit. This is indeed the case. For the integrand of k>, the approximations 
( 3.46D -(3.48) cease to be vahd because we may be of the order one or larger. Treating e as a variable 
and making use of the expansion formula 



cn{-K + e) = cn{K - e) = ^ 1 - — 



V2 



40 1290 



(3.53) 



and the identity 



sn^crdn^cr = — (1 — cn^cr) 
2^ ' 



(3.54) 



we find the approximations of Di{lo), D2{uj) and D±{uj) in terms of I?o(w) and I?i(aj) of the single 
Wilson line, i. e. 



■±ie-K) 



■±ie-K) 



D±(lo) 



\a=±(e-K) 



(3.55) 



and 



a=±(K-e) 



(3.56) 



The correction is of the order which remains small throughout the integration domain of k>. 
The WKB analysis of the appendix A yields 



'ni,2{K - e) = S.i,2{-K + e) = Ci,2(w)ui(a;e) 



1+0 



(3.57) 



and 



il^{K - e) - i+{-K + e) = C3(c^)«i(o.e) 



H 











l + o( - 



(3.58) 



(3.59) 



with C's given by the first two terms of their asymptotic expansions ( 3.27 )-( 3.29 ). Substituting 
eqs.( ^.57 )-( 3.5£ ) into the expression of C/q(cj), we observe that only one term of the numerator of 
(3.4) dominates exponentially. It follows from (3.39) and (3.4C) that 

7^2(c^)=7^?(w)[l + o( -)1 (3.60) 



where we have utilized the relations in ( 3.24 ) and the explicit forms of the functions u's in (3.6)- 
(3.8). Consequently lim^jQ^oo k> = and we arrive at the integral representation of the coefficient 



£1(1) 

3 

7r2 



dw In p{uj) 



(3.61) 



with p{ui) given by (3.50). This integral is well defined and will be evaluated numerically in the 
next section. 
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4. The evaluation of ki 



— - numerical part 



As was explained in section II, the algebraic coordinate z will introduce an artificial singularity 
to the differential equations underlying the determinant ratio of the parallel lines. The conformal 



metric (2.22) we work with involves elliptic functions. This is not a big deal for numerical analysis. 
The elliptic functions can be expressed as the ratios of theta functions [|o| , 



sncT 



and 



where 



^i(z) 



dna — 



oo 



^3^1 








^4^2 


(^) 


•02^4 




^4^3 


im) 


^3^4 


(#) 



2 ^(-)"g("+5)' sin(2n + 1)ttz 

n=0 



oo 



n=0 



cos(2n + l)7rz, 



n=l 



COS 2n,7rz, 



and 



Mz) = l + 2^(-)"g"'cos2n7rz 



(4.1) 
(4.2) 

(4.3) 

(4.4) 
(4.5) 
(4.6) 

(4.7) 



with = i?i(0) for 1 = 1,2, 3, 4. The quantity -ds is related to the complete elliptic integral of the 
first kind via K{k) = ^"d^ with k the modulo. The expansion parameter q — ~ 0.0432139 for our 



k = so the series (4.4)- ([4.7|) converge extremely fast. Also in this case, ??2 



^4 



2-31?.. 



The Schroedinger like equations ( |3.1l| )-(3.13) are solved with the fourth order Runge-Kutta method 
under the boundary conditions ( |3.14|) -( 3.17 ). For eqs.( 3.11 ) and ( [3.12 ), we take advantage of the 



symmetry property (3.18) and evaluate the Wronskian by the formula 

Wl,2i^) = -277i,2(0)77;,2(0) 



(4.8) 



where the prime denotes the derivative with respect to a. The coefficients Ci,2(a;) follows from 
(3.24). For eq.(3.13) with the upper sign, we develop ?7+(o') from a ~ ~K and ^+{cr) from a = K, 



evaluate their Wronskian at ct = and calculate the coefficient 6*3(0;) from eq.(3.24). An alternative 
way is to run the solution 77+ (cr) all the way to K and calculate the Wronskian by eq.(3.25). To 
avoid the rapid changes of the potential function near the singularity a = —K, we start with an 
analytical approximation of 771.2 (cr) and '7+(o') at cr = —K + 5 with 5 « 1 and then run the 
Runge-Kutta iteration for a > —K + 5. Notice that 5 here is not the regularization parameter 
introduced below ( |2.6D and on LHS of ( 3.38| )).On writing x — ujS, we find the approximate solutions 



and 



mM^^ + 6) = 3{ui{x) + ci,2[p{x)ui{x) + q{x)vi{x)]} 



ri+{-K + S) = 3{ui(x) + c+[p{x)ui{x) + q{x)vi{x)]} 



(4.9) 



(4.10) 
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where ui and vi are given by (3.6), 



p{x) 



20a;4 



and 



qix) 



r 



1 



- [ i^x-" + -X + ■ 



— sinh2a;— —x cosh 2a; 



(4.11) 



(4.12) 



The coefRcients ci = 1, C2 = —4 and c+ = — |;. No such a precaution is necessary for the solution 
£,+ {cr) and the Runge-Kutta can start right at the point a = if. The numerical results of Ci{uj), 
6*2(0;) and C3(w) are displayed in Fig.l, where we have introduced 



Ci(c.) = -Ci(c.)e-2^'^ 



and 



-2Kb. 



-2Ku. 



= 1 

LO 

= 1-^ 

2w 



(4.13) 
(4.14) 

(4.15) 



with the last step of each equation following from the asymptotic expansions (3.27), 
(3.29). The comparison of the numerical results with the asymptotic expansions (4.13), 



.28) and 
.14) and 



(4.15) is shown in the table 1 and that with the small w behaviors ( |3.35 ), (3.36) and (3.37) is shown 
in the table 2. The agreement is excellent. To gain more confidence on the numerical solutions 
of the differential equations (3.11), (3.12) and (3.13) for intermediate w, we checked the numerical 
code against a soluble model in which we base the covariant derivatives in ( ^.15 ) on the following 
AdS2 metric 

ds^ = (dr^ + rfcr^) (4.16) 

COS^ cr 

with — f < cr < ^. The differential equations corresponding to ( ^.11 ), ( 3.12| ) and ( 3.13 ) can be 
reduced to hypergeometric equations and the exact forms of the C-coefficients of soluble model are 
derived in the appendix C. They read 



3w^ 



■ sinh 



and 



We have 



and 



4^2 



4^2 + 1 



coshi 



lnp"°'-(cj) EE In 



3iCfi-(cj)4 



Q°i-(w)C|°'-(a;)5 sinh5 ttw 
dwlnp"°'-(w) = 0. 



1 

2^ 



In Fig. 2, we plot the function p(w) of ( 3.50 ) along with that of the soluble model ^''"'■(a; 
small w behavior of the former can be fitted to a polynomial 

piuj) = 2.19171tJ - 3.4445w^ + 6.21735tJ^ - 10.8863a;'^ + 17.5978a;^ 



(4.17) 

(4.18) 
(4.19) 

(4.20) 
The 

(4.21) 



consistent with ( [3.52 ). For large uj, the products uj'^p{lu) and w^psoi. j^^-j g^j.g tabulated in table 3 with 
both determined numerically. Both Fig. 2 and Table 3 suggest that p{uj) falls off faster than p^°^-{uj) 
for large uj. An analytical demonstration requires extending the WKB approximation in appendix 
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LO 


10 


100 


200 


300 


400 


500 


600 


700 




0.848635 


0.847228 


0.847217 


0.847215 


0.847214 


0.847214 


0.847214 


0.847213 


|(i-C2M) 


0.844365 


0.847182 


0.847205 


0.847210 


0.847211 


0.847212 


0.847212 


0.847212 


2c^(l-C3H) 


0.846503 


0.847205 


0.847211 


0.847212 


0.847213 


0.847213 


0.847213 


0.847213 



Table 1: The large u behaviors of the numericaUy generated Ci{ui), C2{io) and C3{uu) 



LO 


0.00001 


0.00005 


0.0001 


0.0005 


0.001 


0.005 


0.01 


Ci(uj) 


10.16655701 


10.16655704 


10.16655711 


10.16655950 


10.16656698 


10.16680621 


10.16755383 




5.08327851 


5.08327852 


5.08327857 


5.08328004 


5.08328465 


5.08343202 


5.08389259 




4.00000000 


4.00000001 


4.00000006 


4.00000143 


4.00000574 


4.00014355 


4.00057424 



Table 2: The small uj behaviors of the numerically generated Ci{u)), Ci{u}) and Ci(tj) 



LO 


10 


100 


200 


300 


400 


500 


600 


700 


Up' lnp(a;) 


-0.000052 


-0.000001 


-0.000002 


-0.000003 


-0.000004 


-0.000004 


-0.000005 


-0.000005 


ij- In iuj) 


0.493798 


0.499938 


0.499984 


0.499993 


0.499996 


0.499997 


0.499998 


0.499999 



Table 3: The large tj behaviors of p(bS) and p°°' (a;). 



A to higher orders and will be rather tedious. Here we merely post our observation without offering 
a rigorous proof. The self-adaptive Simpson integration of In yields 



(iwlnp(ti;) 



-0.56534, 



which, upon substitution into (3.61) leads to 



-1.33460. 



(4.22) 



(4.23) 



The relative error in the numerical valuation of the elliptic functions is about 10^^^ and that of 
the coefficients C^°' (a;)'s extracted from our Runge-Kutta iteration is found below 10^^^. Notice 
that the near singularity expansion of the trigonometric functions pertaining to the soluble model 
goes by the second power of e, while the same type of expansion of the elliptic functions pertaining 



to the parallel lines, ( 3.53 ) and ( 3.54 ), goes by the fourth power of e. Therefore the approximations 
(4.9)-( 4.12 ) should work better for the parallel lines. Likewise is the numerical integration ( 1.22 ), 
the integrand of which vanishes faster than that of the soluble model at large w. For the soluble 
model, we found 2.39x 10~* in contrast to the exact value zero of (4.20). Consequently, the accuracy 



of our numerical algorithm should be amply sufficient for the six significant figures of the k value 
reported in this paper. 



5. Determination of K2 

To determine K2 , we quote two formula of Ref . , one for a bosonic determinant and the other for 
a fermionic determinant. Consider a general 2d metric 

ds^ = gafida°'daf^ (5.1) 
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Figure 1: the top curve represents (73(0;), the middle one represents Ci{ijj), the bottom one represents 





Figure 2: the left curve represents p{uj), while the right one represents p'^°' (cj) 



with the scalar curvature R. Define the functional operator Am = M + X) vi^ith the 



Laplacian with respect to the metric (5.1) and (Af , X) functions of coordinates. Varying M amounts 



to a conformal transformation of the metric (5.1) and associated anomaly contributes a nontrivial 
finite term to the variation of the functional determinant of Am- We have 



In 



detAM) 
det Ai 



1 

An 



InM i-R-X 
6 



1 ^pdlnM d\nM 
12^ dcr^ 



boundary terms, 
(5.2) 

Since we are always taking the difference between the parallel Wilson lines and the two single Wilson 
lines, the boundary terms cancel and we may integrate by part freely. For a Dirac operator with 



respect to the metric (5.1), 7"Vq,, we define 

A^ = -(/C-i7"V„)2+/C-2y 

with /C and Y functions of coordinates. The measure transformation formula corresponding to ( |5.2D 
reads. 



(5.3) 



^ detA£)\ 1.2^ 



ln/C(-i? + 2r +-5"^^^ 



boundary terms, (5.4) 



where we have multiplied the integral in [|l^ by two, taking into account that A^ here is a 2 x 2 
matrix in the spinor space. 



Coming to the determinants we are interested in, metric (2.25)and metric (2.22) are all confor- 
mal with 

gai3 = e^^^Sap (5.5) 
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and the scalar curvature 



We have (ctq, <ti) = (r, and x = InC for the single line, and (croiCn) = (t, cr) and x = ln(-\/2cno-) 
for the parallel lines. With the measure scaling functions M — e^^ and /C = e-^, those functional 



operators of (2.27)-(2.33) without hats corresponds to Ai and Af of eqs.(5.2) and ( ^.4[ ) and that 
with hats to Am and A^ there. The "mass square" X of (5_^) equals to zero for Ao[Ci] and A0P2], 



equals to 2 for Ai[Ci] and A 1^2] and equals to 4 + i? for A2[C2]. The "mass" Y of (^J) equals to 
one for all fermionic determinants. It follows from (5.2) that 



In 



In 



det Ao[Ci] 
det Ao[Ci]^ 

det A0P2] ' 
det A0P2] 



fin. 



fin. 



1 

12^ 



aa dx dx 



boundary terms, 



In 



In 



det Ai[Ci] 
det Ai[Ci]^ 

detAi[C2]^ 
det Ai[C2] ) 



= ^ / d^f^S"^ S^S^ + boundary terms, 

d ae '^x + boundary terms, (5.9) 



(5.7) 
(5.8) 



127 



and 



In 



detA2[C2] 
det A2[C2] 





p. 9x 


dx 






J dx 


dx 



'C2 



^^ae + boundary terms. (5.11) 



fin. 



'C2 



where the subscript of the integration sign indicates the world sheet integration extends to. Simi- 
larly, the formula (^.4[) implies that 



In 



and 



In 



|det£iF[Ci 
I det Df[Ci 



I det Df[C2 
\detDF[C2 



fin. 



A_ f d'^s-^^p,^!- [ d 



'^cre ^-^x + boundary terms (5.12) 



1 

24^ 



a/3 9x dx 



— / d^ae ^-"^x + boundary terms. (5.13) 



Substituting into ( 2.35 ) and ( ^.37| ) for the single line and the parallel lines, we find their contribu- 
tions add up to zero in each case i.e. W2P1] — W2P2] — 0. Consequently, 



K2 = 0. 



(5.14) 



This, together with (4.23) leads to our final result (h 



6. Concluding remarks 

As AdS/CFT has become an important reference to understand the observation of the strongly 
interacting quark-gluon plasma created by heavy ion collisions, it is critical to asses the robustness 
of the leading order prediction by exploring the next order correction in the expansion according 
to the inverse powers of the large 't Hooft coupling A — A^c<?ym- "^bc subleading terms of the 
expansion have been addressed in the literature in the context of the equation of state and the 
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shear viscosity . This type of corrections comes from the a'^ correction of the target space 

metric Its contribution is of the order 0(A~'^/^) relative to the leading order in the Af — 4 

SYM and is present only at nonzero temperature. In case of the expectation value of a Wilson loop 
operator, however, the dominant correction stems from the fluctuation of the world sheet around 
its minimum area and is suppressed only by 0{X~^^'^) relative to the leading order. It shows up at 
all temperatures and is more difficult to compute. The only attempts made in the literature in this 
regard include the strong coupling expansion of a single line, a circular loop and a spinning line at 
zero temperature [|l^[^[0[0- These Wilson loops, though theoretically important, do not carry 
direct phenomenological implications. 

In this work, we have extended the method in [|l^ to the fluctuations of the world sheet dual to 
a pair of parallel Wilson lines and have derived the next term of the strong coupling expansion of the 
heavy quark- ant iquark potential inAf — 4 SYM at zero temperature. We start with the determinant 
ratio for a single Wilson line and that for parallel lines in the static gauge, in which the fluctuations 
come from eight transverse bosonic coordinates and eight 2d Majorana fermions. Then we scaled 
the operators underlying the determinants, leaving a trivial measure for the associated spectral 
problem. The subleading term of the heavy quark potential is extracted from the combination 



(2.17), which consists of the spectral and the measure parts. A robust numerical result of the 
former is obtained and the contributions from measure change of each determinant cancel. We 
have, 

«(A) 



V{r) 



with 



47r2 V\ 



a(A) = 



A 



1 - 



1.33460 



1 - 



27r2 



In- 



Va 

27r 



for A >> 1 



7B + 1 +0(A') 



(6.1) 



(6.2) 



for A << 1 



where the weak coupling expansion obtained in [|l3[ |29|| from held theory is also included for com- 
pleteness. The authors of also worked out the strong coupling expansion under the ladder 
approximation in field theory. 



ladder 



(r) 



Va 



1 - 



TT 

7t 



(6.3) 



It is interesting to notice that our subleading term is of the same sign as theirs but the magnitude 
relative to the leading order is smaller in our result. In view of the range of the 't Hooft coupling 
which was used for the RHIC phenomenology. 



5.5 < A < 6tt 



(6.4) 



the correction to the leading order of the strong coupling may be significant in magnitude. One 
may define an effective coupling 

VA' = %/A- 1.33460 (6.5) 



If A of (3.4) is replaced by A', the range of the 't Hooft coupling is shifted to 

13.54 < A < 32.22 



(6.6) 



At a nonzero temperature T, however, the order 0(A ^Z^) is not merely a redefinition of the coupling 
and the strong coupling expansion of the heavy quark potential becomes 



V{r) 



4^2 7a 



goirT) - 



1.334605i(rT) 

7x 



(6.7) 
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ds^ = — fiz)dt^ + dx'^ + -j^dz^ + dnl (6.8) 



with 50(2;) and gi{x) two functions satisfying the conditions go(0) — gi{0) — 1. The function .go (a;) 
have been determined by the minimum area of the world sheet in the Schwarzschild-AdSs x target 
space [0 

1 ^ . . ... J_ 

with f{z) = 1 — TT^r^z"* and t the EucUdean time. The one loop effective action underlying the 
function gi (x) has been developed in and the methodology employed in this work can be readily 
generalized there. 

While simple in practice, the static gauge we worked with suffers a problem. Though the 



combination (2.17) gives rise to a finite result, neither the UV divergence nor the conformal anomaly 



of each term on RHS of ( 2.17 ) vanishes. A less problematic gauge is the conformal gauge, in which 
the world sheet metric is not set to the induced metric at the beginning. One has to include the 
determinant of the longitudinal fluctuations and that of the ghost and an appropriate measure of 
the path integral. The contributions from the transverse bosons and fermions obtained in this paper 
will remain there, but other contributions including the measure change may be subtle to collect. 
It is important to carry out the parallel analysis in the conformal gauge to ascertain that our result 
in this paper is complete. Another alternative is the canonical quantization method employed in 
II31II . We hope to report our progress in this direction in near future. 
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A. 

To extract the large w behavior of the coefficients Ci{uj), C2(w) and C3(w), we introduce x = 
uj{K + a) and y = uj{K — a). For a + K « 1{K — a « 1), the solutions of the differential 
equations can be approximated by that of the equations for the single line, which extends to 
X >> l{y >> 1) for large w. The WKB approximation applies for a; >> 1 and y >> 1. In case 
of eq.( p.l3| ) with upper(lowcr) sign, the WKB solution can be extended all the way to the point 
a — K{<7 ~ —K) and the requirement y » \{x » 1) may be relaxed. We match the single line 
solution and the WKB ones in the regions where both approximations apply. 



Consider the equation (3.11) first. We start with the approximate solution near a = —K 

771 ~ 3^1(2;) (A.l) 



with a + K << 1, where the coefficient 3 follows from the requirement ( 3.14 ). The asymptotic form 
for x >> 1 reads 

The WKB solution to be matched is given by 



7;i~exp / da'\/c^2+ . (A.3) 
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Expanding the square root for large to and using the derivative formula 

d sncrdncr 1 / 1 



da 



2 V cn^cr 



cn^cr 



we find that 



, 1 

?7i ~ /I exp Loa + — 



da' 



sncrdncr f ; , / 2 / 
A uja H — / da en a 



wcncr 2uj 



K 



(A.4) 



(A.5) 



with A a constant to be determined. In the left matching region where x >> 1 and a + K « 1, 
the approximations 

sn(cr)dn(cr) 1 



cn(cr) 



K + a 



(A.6) 



and dcr'cn^cr' ~ yield the coefficient A = |e^'^. In the right matching region where K — a« 
1 and y >> 1, the WKB solution (A^) becomes 



3 2Kuj--^ -v+- 

r]i ~ 776 e , 



(A.7) 



where we have used the approximation 



sn(cr)dn(CT) 
cn(cr) 



K-a 



there and the constant 



1 

'^2 



K 



dacv?a ~ \f2 \ dx 



-K 



Comparing with the expression of (3.6), we obtain that 

3 



(A.i 



(A.9) 



(A.IO) 



The asymptotic behavior of ( ^.27 ) is extracted in the limit y — >■ and the relation ( 3.57 ) for 771 and 
^1 follows. 

The equation (3.12) can be treated similarly. We start with the same expression of (A.I) for 
7/2(0-) near a = —K but replace the WKB solution (A. 3) by 



7?2=^exp^y da'^juj^ + -^-cv?a^ 



(A.ll) 



Eqs.(AJ) and (A_/7) become 



m ^ A exp [ uja + — 



with the same A and 



K 



da' 



cn^cr' 2uj 



—— I da'cn^a' 



K 



A exp Loa 



sncrdncr 1 



J / 2 / 
da cn a 



K 



I \ '3 2Ku-^ -v+- 



(A.12) 



(A.13) 



for y >> 1 and K — a « 1. The asymptotic behavior ( 3.28 ) and the relation ( 3.57 ) for (772,^2) 
are extracted then. 



Coming to eq.(3.13), the single line solution ( [A.lj ) remains approximating and we have 

ry+~3ui(a;) (A.14) 
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for a + K << 1. The WKB solution it matches with for 2: >> 1 reads 



77_|_ ~ exp / da' i ui'^ 



sna'dncr' 



3 / ^ 1 r , , 1 - VZsno-'dna' 
- exp Kuj + uja + — da — 



— exp Koj + uja ^ 
2 V 2v^a; cna 



1 v^sncidncr — 1 ^ j / 2 / 



da cn a 



-K 



and works all the way to the point a = K . Near that point, we find 

3 



2Ku 



(A.15) 



(A.16) 



Eqs.(|3.29p and ( |3.58| ) follow from the form of wo(y), ( |3.24| ) and ( ^.25^ . The relation ( |3.59D is obtained 
starting with the WKB solution 



sinh 



K 



da\ Lo"^ + 



sncr'dncr' 



2cn2cr' V^cn2cr' j 
and matching it to the approximate solution vi{x) for K + a « 1 and x >> 1. 



(A.17) 



B. 



The differential equation (3.11) at uj = can be converted to a hypergeometric equation by the 
transformation 



X = cn a 



[a) = Va;/(a;), 



I.e. 



We have 



ryi(cr) = 2cn2o-i^ ( i, ^; ^;cn''cr 



and ^i(a) = 771 (—fi). It follows from the formula 



^F{a, h; c; z) = —F{a + 1, 6 + 1; c + 1; z), 
dz c 



F(a,6;c;l) 



T{c)T{c - a~b) 
r{c~a)T{c-b) 



for Re(c — a — fe) > and 



r(a)r(6) 

for Re(c ~ a ~ b) < and e > that the Wronskian 

48 r(|)r(i)r(^)r(i) 



W[m,^i] = -2 hm ^i(-e)^;(-e) 
727ri 



7v/2r(|)r(i) r(|)r(i) 



(B.l) 
(B.2) 

(B.3) 

(B.4) 
(B.5) 

(B.6) 



(B.7) 
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Divided by -3, we derive ( ^.35 ). 



With the same transformation (B.l), cq.(3.12) becomes 

d'^f (7 9 \ df 1 



We have 



4 A J dx 4 



1 7 



r]2{a) = 2cn aF 1, -; -;cn cr 



4' 4' 



(B. 



(B.9) 



and ^2(0") = 772(— cr). It follows from (B.4)-( |B.6| ) that the Wronskian 

32 r(|)r(i)r(^)r(i] 



iyfe,6] = -2 U-eW,{-e) = r (I) r (i) r(2)r (|) 

367r5 



(B.IO) 



and ( 3.36 ) is obtained as — W[?72, C2]/3. Using the series representation of the hypergeometric 
function, we find that 



771 — —X 4 
' 2 



dx'x'-^ _ 3%/2 



dcr'cn^cr'. 



(B.ll) 



if 



But we fail to find a similar expression for 772(c). 

As to (73(0), we notice that the solution of the 1st order differential equation 

V2cna^ + = 
da 



(B.12) 



also solves eq.( |3.13 ) with the upper sign. The eq.(|l|) can be solved readily and we obtain 

V'(cr) = B 



'V2d 



na — sncT 



V2dncr + sncr 

with B a constant, where we have used the indefinite integral 



1 1 \/2dna — sncr 
da = = In — — h const. 



cncr y/2 V2dncr + sncr 



(B.13) 



(B.14) 



as can be verified by taking derivatives of both sides. Setting the constant B = 2, we find that the 
function ip{a) satisfies the boundary condition of f+(cr) at a = K and therefore (cr) = '0(o')- As 

cr ^ -K, 

4 



a + K 



(B.15) 



and we end up with C3(0) = 4. 



In this appendix, we present the details of the soluble model which is introduced to check our 
numerical algorithm. The model is largely motivated by the work in We shall use the same 
symbols (ry,^) for the solutions of the counterparts of the differential equations ( ^.ll| ), ( [3.12 ) and 
( 3.13|) . Because the scalar curvature of the metric ( 4.16 ) is R = —2, the counterparts of ( 3.11 
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and ( 3.12 ) are the same. Consequently, -Di(aj) — D2{i^) and (»7i,Ci) = ('72,^2) in this case. The 
counterpart of the eq.( 3.11 ) or ( 3.12 ) reads, 

(P6 



(C.l) 



and has the same set of indexes at the regular points a — ±^ as that of ( ^.11 ). The symmetry 
property (3.18), the relation ( 5.24) and the formula (4.S) remains valid. The solution ?7i((t), specified 
by the boundary conditions ( 3.14 ) with K replaced by ^ is 



771 (cr) = (wcos cr)^F + ^, ^, 



(C.2) 



and ^(cr) = 7y(— cr). It follows from (4^), the formula (B.4) 

a;4(a;2 + 4) 



~2 lim ?7i(-e)77i(-e) 

£->0+ 



for hypergeometric functions that 

r(|)r(|)r2(i) 



5 r(2±-)r(2 + zf)r(i^)r(2-zf 



(C.3) 



1 



■ sinh vroj 



(C.4) 



Divided by —3a; we end up with (4.17). The nonzero component of the spin connection corre- 
sponding to the metric ( 1.16 ) is = tana and the counterpart of eq( |3.13 ) with the upper sign 
reads 

(P(h f ^ 1 



1 + sin cr 



= 0. 



(C.5) 



The equation (C.5) can be reduced to a hypergeometric equation and the solutions satisfying the 
boundary conditions (B.14), (3.16) and (3.17) (with K replaced by ^) read 



?/-)- = 2uj (1 + sin cr)i<" 1 + icj, 1 — iuj 



5 1 + sin cr 



and 



i+ = -y|(l + sincr)\/l - sin crF ( - + - - iw; -; 



3 1 — sin cr 



(C.6) 
(C.7) 



Their Wronskian 



r(l)r(i] 



3^2 



The eq.(p^) follows then from (OJ) and (|^). 



cosh? 



(C.8) 
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